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Abstract — In this paper, a new filtering method to remove 
Rician noise from magnetic resonance image is presented. This 
filter is based on Neutrosophic set (NS) approach of median 
filtering. A Neutrosophic set, a part of neutrosophy theory, 
studies the origin, nature and scope of neutralities, as well as 
their interactions with different ideational spectra. Here the MRI 
image is transformed into NS domain, which is described using 
three membership sets: T, I, F. The entropy of the neutrosophic 
set is defined and employed to evaluate the indeterminacy. The y- 
median filtering operation is used on T and F to decrease the set 
indeterminacy and remove noise. We have conducted 
experiments on real MR image with Rician noise added. The 
performance of this filter is compared with the median filtering 
and the classical non local mean (NLM) de-noising approaches. 
This filter outperforms the median filter for different levels of 
noise. For low signal to noise ratio (SNR), this filter provides high 
peak signal to noise ratio (PSNR) than the NLM approach. 

Keywords- Denoising, entropy, magnetic resonance imaging, 

neutrosophic set, PSNR, Rician distribution . 



I. Introduction 

Magnetic Resonance Imaging (MRI) is a powerful 
diagnostic technique used in radiology to visualize detailed 
internal structures of the human body. MRI is based on the 
principles of nuclear magnetic resonance (NMR), to map the 
spatial location and associated properties of specific nuclei or 
protons in a subject using the interaction between an 
electromagnetic field and nuclear spin [1]. MRI is primarily 
used to demonstrate pathological or other physiological 
alterations of living tissues and is a commonly used form of 
medical imaging. Despite significant improvements in recent 
years, MR images often suffer from low signal to noise ratio 
(SNR) especially in cardiac and brain imaging. This is 
problematic for further tasks such as segmentation of 
important features; classification of images for computer aided 
diagnostics, three dimensional image reconstruction and 
image registration. Therefore, noise reduction techniques are 
of great interest in MR imaging. 



The MR image is commonly reconstructed by computing 
the inverse Discrete Fourier Transform of the raw data. The 
signal component of the measurement is present in both real 
and imaginary channels. Each of the orthogonal channels is 
affected by additive white Gaussian noise. The noise in the 
reconstructed complex valued data is thus complex white 
Gaussian noise. Most commonly, the magnitude of the 
reconstructed MRI image is used for visual inspection and for 
automatic computer analysis. Since the magnitude of the MRI 
signal is the square root of the sum of the squares of two 
independent Gaussian variables, it follows Rician distribution 
[2]. The Rician noise is signal dependent and is therefore 
difficult to separate from the signal. 

Numerous approaches of denoising MR images have been 
proposed including approaches based on anisotropic diffusion 
[3-7], the wavelet transform [8-13], bilateral and trilateral 
filtering [14-16], the non-local means algorithm [17- 22] in the 
literatures. In this paper, the neutrosophic set approach (NS) 
filter for image denoising proposed in [23] is adapted to 
remove the Rician noise in MRI. First the noisy MRI is 
transformed into the neutrosophic set and the y -median 
filtering is employed to reduce the indetermination degree of 
the image, which is evaluated by the entropy of the 
indeterminate subset, after filtering, the noise will be removed. 
The experiments on real MR images having different Rician 
noise levels demonstrate that this approach can perform 
denoising better. 

II. Neutrosophic MR Image Denoising 

Neutrosophy, a branch of philosophy introduced in [24] as a 
generalization of dialectics, studies the origin, nature and 
scope of neutralities, as well as their interactions with different 
ideational spectra. Neutrosophy theory considers proposition, 
theory, event, concept or entity, <A> is in relation to its 
opposite <Anti-A> and the <Neut-A> which is neither <A> 
nor <Anti-A>. The neutrosophy is the basis of the 
neutrosophic logic, neutrosophic probability, and neutrosophic 
set and neutrosophic statistics [24]. In neutrosophic set, the 
indeterminacy is quantified explicitly and the truth- 
membership, indeterminacy-membership and falsity- 
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membership are independent. The neutrosophic set is a general 
formal frame work which generalizes the concept of the 
classic set, fuzzy set, interval valued fuzzy set, intuitionistic 
fuzzy set, and interval valued intuitionistic fuzzy set, 
paraconsistent set, dialetheist set, paradoxist set and 
tautological set [24]. The neutrosophic set had been applied 
into image processing such as denoising [23], segmentation 
[25, 26]. The definition of a neutrosophic set and its properties 
are described briefly. 



S(i, j ) = abs(g(i, j ) - g(i, j )) 


(4) 


F(i,j) = \-T(i,j ) 


(5) 



where g(i,j) is the local mean value of the pixels of the 
window. S(i, j) is the absolute value of difference between 
intensity g(i, j ) and its local mean value g(i , j) . 



A. Neutrophic Set 

Definition 1 (Neutrosophic Set). Let U be a Universe of 
discourse and a neutrosophic set A is included in U . An 
element x in set A is noted as x{T,I,F). T,I,F are real 

standard and non standard sets of ] - 0, 1 + [ with 
sup T = /_sup, inf T = t _ inf, sup/ = z sup, inf / = / inf, 
supF = /_sup,inf F = f _ inf and ?/_sup = /_sup + 
/_sup + /_sup, n _mf = /_inf +/_inf+ / inf. 

T, I and F are called the neutrosophic components. The element 
x(r,/,F) belongs to A in the following way. It is t% true in 
the set, /% indeterminate in the set, and f% false in the set, 
where t varies in T , / varies in / and / varies in F . 

B. Transform the image into neutrosophic set 

Definition 2 (Neutrosophic image). Let U be a Universe 
of discourse and IF is a set of U , which is composed by 
bright pixels. A neutrosophic image P NS is characterized by 
three membership sets T,I,F. a pixel P in the image is 

described as P(T, I, F) and belongs to W in the following 
way: It is t true in the set, i indeterminate in the set, and 
/ false in the set, where t varies in T , i varies in / and 
/ varies in F . Then the pixel P(i , j) in the image domain is 
transformed into the neutrosophic set 
domain P NS ( i,j) = {T(i,j),I(i,j),F(i,j )} . T(i,j),I(i,j ) 
and F(i,j) are the probabilities belong to white pixels set, 
indeterminate set and non white pixels set respectively, which 
are defined as: 

T(i, J) = MUhjgjn. (1) 

S max — S min 



C. Neutrosophic image entropy 

For a gray image, the entropy is utilized to evaluate the 
distribution of the gray levels. If the entropy is the maximum, 
the intensities have equal probability. If the entropy is small, 
the intensity distribution is non-uniform. 

Definition 3 (Neutrosophic image entropy). Neutrosophic 
entropy of an image is defined as the summation of the 
entropies of three subsets T, I and F : 



En NS = En T + Enj + En F 


(6) 


maxjT 7 } 

En T = - ^p T (i) In p T (i) 


(7) 


i=min{T} 

max{/} 

Enj=- y Pl (i) In p/(i) 


(8) 


/=min{/} 

max{F} 

En F =- y Ph (i) In p F (i) 


(9) 



/=min{F} 



where En T , Enj and En F are the entropies of sets T , / and 
F respectively. Pt(i) > Pi ( 0 and Pp{i) are the 
probabilities of elements in T 7 ,/ and F respectively, whose 
values equal to i . 

D. J - median filtering operation 

The values of /(/, j) is employed to measure the 
indeterminate degree of element/^ ifj) • To make the set 

/ correlated with T and F , the changes in T and F 
influence the distribution of element in / and vary the entropy 
of / . 

Definition 4 ( y - median filtering operation). A y — 

median filtering operation for P NS ’Pns(y) , is defined as: 



i i+w/ 2 

s(U)= I 

WXW /o 

vv vv m-i-w/2 


j+w/2 

Z g( m > n ) 

n=j-w/2 


(2) 


Es (r) = P(T(y), i{r), F(j)) 






(3) 


f( r )=l T 1<r 

’ \Ty I>r 



( 10 ) 

(ii) 
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T r {i,j)= median {T(m,n)} 

( m,ny=Sij 


(12) 


mA F - ,<r 

’ \F r i>r 


(13) 


F r (/, /)= median {F(m, n)} 

(m,n)eS i j 


(14) 


lJi,j)= 

rV J> Sf -Sf . 

T max T mm 


(15) 


S f ( i , j) = abs(T(i , j) - T (; i , j )) 


(16) 


— i i+w/2 j+w/2 

T{i,j)= X X T(m,n) 

\/ nil 


(17) 



wXw m=i-w/ 2 n= j-w/ 2 



where Sj. (/, j) is the absolute value of difference between 

intensity f(i,j) and its local mean value T(i,j) at 
(i, j) after y - median filtering operation. 

The new neutrosophic approach to image denoising is 
described as below: 




Step 1: Transform the image into NS domain; 

Step 2: Use y - median filtering operation on the true subset 

T to obtain Ty ; 

Step 3: Compute the entropy of the indeterminate subset 

Iy , En * (i) ; 

En f (i + 1 )- En- (i) 

Step 4: if — < 8 , go to Step 5; 

En f (7) 

Else T = fy, go to Step 2; 

Step 5: Transform subset Ty from the neutrosophic domain 
into the gray level domain. 

III. Results and Discussion 

In the experiments, we used the T2 Axial with TR 9000, 
TE 87 of 5th Slice and T1 Sagittal with TR 552, TE 17 of 5th 
Slice of brain MR image of 2 years old male child, which is 
obtained from SIEMENS 1.5 T MRI scanner. The simulated 
Rician noise with different level is added to the both the 
images. 

The performance of the denoising algorithm is measured 
by using the peak- signal-to -noise ratio (PSNR). The higher the 
PSNR is, the better the denoising algorithm is. 

First, the NS filter is compared with the median filter in 
removing the Rician noise with different noise level from the 
both T2 Axial and T1 sagittal brain MR images. Then, the NS 
Filter is compared with the Non Local Mean (NLM) filter 
proposed by Buades et al [17]. 



(c) (d) 




(e) 

Fig. 1 a) Original Image (T2 Axial Brain MRI), b) Noisy image with PSNR 
24.79, c) Denoised using Median Filter PSNR 25.96, d) Denoised using 
NLM PSNR 28.36, e) Denoised using NS Filter PSNR 26.86 

For T2 Axial brain MRI, Fig. 1 (a) is the original image of 
the T2 Axial brain MRI, and Fig. 1 (b) is the image resulted 
from the Fig. 1 (a) added with Rician noise (s=20; 
PSNR=24.79). Fig. 1 (c), Fig. 1 (d) and Fig. 1 (e) are the 
denoised results of Median filter (PSNR=25.96), NLM 
(PSNR=28.36) and NS Filter (PSNR=26.86) respectively. Fig. 
2 shows the PSNR comparison of Median and NS filter. From 
these results, we conclude that the NS filter performs better 
than the Median filter for different noise level. Fig. 3 (a) and 
(b) shows the PSNR comparison of NLM and NS filter. From 
these performance results, we can conclude that at High SNR 
(low noise level), NLM performs well than NS filter and at 
low SNR (high noise level), NS filter performs better than that 
of NLM. 
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Median Vs NS Filter 




Fig. 2 The performance comparison of Median and NS Filter on T2 Axial 
brain MR Image with Rician noise. 



NS Filter Vs NLM: High SNR 




Noise 



(a) 



NS Filter Vs NLM : Low SNR 




Noise 

(b) 

Fig. 3 The performance comparison of different filters on T2 Axial brain MR 
Image with Rician noise, a) High SNR, b) Low SNR. 



For T1 Sagittal brain MRI, Fig. 4 (a) is the original image 
of the T1 Sagittal brain MRI, and Fig. 4 (b) is the image 
resulted from the Fig. 4 (a) added with Rician noise (s=20; 
PSNR=26.80). Fig. 4 (c), Fig. 4 (d) and Fig. 4 (e) are the 
denoised results of Median filter (PSNR=27.19), NLM 
(PSNR=29.22) and NS Filter (PSNR=28.19) respectively. Fig. 
5 shows the PSNR comparison of Median and NS filter. From 
these results, we conclude that the NS filter performs better 
than the Median filter for different noise level. Fig. 6 (a) and 
(b) shows the PSNR comparison of NLM and NS filter. 




(c) 



(d) 




(e) 

Fig. 4 a) Original Image (T1 Sagittal Brain MRI), b) Noisy image with PSNR 
26.80, c) Denoised using Median Filter PSNR 27.10, d) Denoised using 
NLM PSNR 29.22 e) Denoised using NS Filter PSNR 28.23 
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Median Vs NS Filter 




Fig. 5 The performance comparison of Median and NS Filter on MR Image 
T1 Sagittal with Rician noise. 

NS Filter Vs NLM : High SNR 




Noise 

(a) 



NS Filter Vs NLM : Low SNR 




(b) 

Fig. 6 The performance comparison of different filters on T1 Sagittal MR 
Image with Rician noise, a) High SNR, b) Low SNR. 



From these performance results, we can conclude that at High 
SNR (low noise level), NLM performs well than NS filter and 
at low SNR (high noise level), NS filter performs even better 
than that of NLM. 

From Table I, we better understand that the NS filter 
perform well than Median filter for different noise level and 
also, for low SNR the NS filter performs better than NLM 
irrespective of the MR images. 



TABLE I 

C OMP ARISION OF MEDIAN, NLM, NS FILTER BASED ON PSNR IN dB 





PSNR (dB) ( 


Noise 

Level 


Median 


NLM 


NS filter 




T2 

Axial 


T1 

Sagittal 


T2 

Axial 


T1 

Sagittal 


T2 

Axial 


T1 

Sagittal 


2 


23.18 


27.75 


35.65 


46.56 


32.69 


33.40 


10 


23.14 


27.24 


34.36 


35.23 


30.41 


31.42 


20 


22.79 


25.89 


28.36 


29.19 


26.84 


28.14 


30 


22.04 


24.20 


24.81 


25.58 


24.19 


25.59 


40 


21.09 


22.48 


21.30 


22.90 


21.99 


23.28 


50 


19.94 


20.86 


20.27 


20.80 


21.05 


22.00 


60 


18.81 


19.44 


19.64 


19.07 


20.02 


20.36 


70 


17.74 


18.08 


18.25 


17.55 


19.90 


19.11 



IV. Conclusions 

This paper presents a neutrosophic denoising approach 
based on neutrosophic set for removing Rician noise from 
brain MRI. This filtering method tends to produce good 
denoised image not only in terms of visual perception but also 
in terms of the quality metric (PSNR). This filter performs 
better than Median filtering method for removing the Rician 
noise with different noise levels. Further, and also it 
outperforms the Non Local Mean approach when the noise 
level is high (low SNR). 
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